Novel method for direct dosimetry of radioactive nuclei in different types of samples using solid-state detectors of nuclear traces

ABSTRACT

The invention concerns a novel method adapted to experimental conditions, for determining the function defining the actual detecting sensitivity of a solid-state detector of nuclear traces and which represents the detecting performances for each type of nucleus emitting energy alpha particles Ei, located in each position inside a sample to be analyzed, and hence precision dosimetry of the radioactive nuclei in a sample using nuclear trace solid-state detectors. The method provides an improvement of the quality of sample analyses, a significant reduction of analysis costs and an important widening of the scope of application of this type of detectors. Thus, the method will bring about the development of a new generation of detectors having a better price-quality ratio which can be found in the field of radioactive nuclei analyses in different types of samples.

RELATED U.S. APPLICATIONS

Not applicable.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

Not applicable.

REFERENCE TO MICROFICHE APPENDIX

Not applicable.

FIELD OF THE INVENTION

The present invention relates to a method for direct dosimetry of radioactive nuclei in different types of samples using solid-state detectors of nuclear traces.

BACKGROUND OF THE INVENTION

The history of solid state nuclear track detectors began in 1958 after initial observations by D. A. Young of fission fragment tracks of lithium fluoride (LiF) irradiated by thermal neutrons. In 1959, Silk and Barnes made these same observations on mica using an electron microscope. Numerous observations have been made since then on a large number of materials. The detectors have the ability to record the passage of particles such as alpha particles, protons, and fission fragments.

In the years that have passed, work by the international scientific community on solid state nuclear track detectors made it possible to develop several dosimetry methods for radioactive nuclei with the use of these detectors. Though these methods use basic knowledge of calculation of the probability that an alpha particle emitted from a radioactive nucleus of a sample can reach and be recorded on a solid state nuclear track detector, they are based on estimation of the global critical angle depending on the sample, the characteristics of the particles to be detected, and the detector. The estimation thus induces an error in the calculation and makes the solid state nuclear track detectors (SSNTDs) imprecise.

Other methods of detection that do not use solid state nuclear track detectors already exist. They involve, for example, geochemical analyses and instrumental techniques such as analysis by neutron activation, analysis by mass spectometry or analysis by gamma spectrometry and liquid scintillation analysis. All of these techniques, often destructive to the samples analyzed, use equipment that is heavy, sophisticated, and especially costly.

BRIEF SUMMARY OF THE INVENTION

In order to correct these disadvantages, a new method has been developed. It is based on the implementation of a new technique specified for calculating the sensitivity of a solid state nuclear track detector for each type of nucleus emitting alpha particles having energy E₁ located in each position in the sample to be analyzed. This method does not allow the estimation but instead the calculation without errors of the critical angle of incidence depending on the sample, the characteristics of the alpha particles to be detected, and the detector. Thus, this method is based on the actual function of the critical angle of incidence which makes it possible to solve the problem of the sample-detector geometry, and as a result, determination of the function that defines the real sensitivity of detection of a solid state nuclear track detector (SSTND) and which represents the efficiency of detection for each type of nucleus emitting alpha particles having energy E₁ located in each postion in the middle of the sample to be analyzed.

Applications of this method allow an improvement in the quality of sample analyses, a considerable reduction of costs for analysis and a great enlargement in the field of application for this type of detector. Thus, a new generation of detectors will be born having a better price-to-quality ratio than has existed until now in the domain of analyses of radioactive nuclei in different types of samples.

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

FIG. 1

FIG. 2

FIG. 3

FIG. 4

FIG. 5

DETAILED DESCRIPTION OF THE INVENTION

II. Methodology and Procedures for Evaluation of Different Parameters

Let (dv) be the volume element of a cylinder having a radius between r and r+dr and a height between h and h+dh inside the sample analyzed, see FIG. II-1 in the drawings section where: (1: SSNTD, 2: plastic capsule, 3: sample analyzed). dv=2πrdrdh

Let dN_(i) be the number of the alpha particle having energy E_(i) emitted by the radioactive element having index i contained in the volume dv and for the exposure time te: dN _(i) =A _(ci) ^(te) ·dv

Where A_(ci) ^(te) is the alpha activity of the radioactive nuclei of the emitter element i of the alpha particle having energy E_(i) per unit volume and for the exposure time te.

It is thus necessary to multiply the number dN_(i) of alpha particles emitted in the volume dv by a probability of detection P_(i)(r,h), [P_(i)(r,h) is the probability that an alpha particle of energy E_(i) emitted by the nucleus of index i contained in the volume dv can be detected by the solid state nuclear track detector] and which also represents the efficiency of detection of alpha particles of energy E_(i) emitted by the nucleus of index i contained in the volume dv.

As a result, for radioactive nuclei of element i contained in the volume dv, the number of tracks developed will be: dn _(i) =A _(ci) ^(te) ·P _(i)(r,h)·dv

Where: dn_(i) is the number of alpha particles having energy E_(i) emitted by the nuclei of index i contained in the volume dv, for the time of exposure (te) and which are detected in our solid state nuclear track detector.

Actually, the total number of tracks observed due to the radioactive nuclei of index i contained in the total volume is: n _(it)=2π·A _(ci) ^(ie) ·∫∫P _(i)(r,h)·rdrdh II-1 Method for Calculating the Probability P_(i)(r,h).

Let a nucleus of index i denoted K_(i), an emitter of alpha particles of energy E_(i), be located in a position of coordinates (r,h) (see FIG. II-5 in the drawings section) such that these emitted particles give rise to observable tracks on the solid state nuclear track detector. The particles must satisfy two conditions:

-   -   1—they are emitted at a well-defined angle of incidence by means         of the function of the critical angle of incidence (θ_(ic)).     -   2—they reach the detector.         II-1-1 Identification and Evaluation of the Function (θ_(ic)).

A) Identification of the Function (θ_(ic)).

FIG. II-2 in the drawings section shows the experimental and intrinsic function of a SSNTD describing the variation of the critical angle (θ_(c)) resolving the tracks as a function of the residual energy for conditions enabling chemical resolution and reading of well-defined tracks. This function is from a solid state nuclear track detector that is currently the most used (CR-39; LR-115). Three possible cases can be distinguished: *E_(ind)<E_(min), *E_(ind)≦E_(min)≦E_(min), *E_(ind)>E_(max) where:

-   -   E_(ind) is the residual energy at normal incidence on the         detector of the alpha particle having energy E_(i) emitted from         a nucleus located at a height h from the plane of the detector.     -   E_(min) is the minimum energy detectable at normal incidence.     -   E_(max) is the maximum energy detectable at normal incidence.

-   Comment: The identification of (θ_(ic))=f(h), (θ_(ic1))=f_(a)(h),     and (θ_(ic2))=f_(b)(h) can also be made by means of the inverse     function of the function (θ_(c))=f(E_(res)):E_(res)=f(θ_(c)) which     shows the minimum and maximum energies detectable for each incidence     (θ_(i)).     A-1—Case of Nuclei Emitters of Alpha Particles Having Energy E_(i)     and with E_(ind)<E_(min),

These are nuclei located in the volume for contribution to the tracks defined by h∈[R_(i)−R_(min),R_(i)] for the case without an absorbent medium, given R_(i) and the total trajectory in the sample of the alpha particle having energy E_(i) in the sample, and given R_(min) and the total trajectory in the sample of an alpha particle having energy E=E_(min). All alpha particles having energy E_(i) emitted from these nuclei will reach the detector with an energy lower than E_(min) and as a result will not be detected.

A-2—Case of Nuclei Emitters of Alpha Particles Having Energy E_(i) and E_(min)≦E_(ind)≦E_(max).

These are nuclei located in the volume for contribution to the tracks defined by h∈[R_(i)−R_(max), R_(i)−R_(min)] for the case without an absorbent medium, given R_(max) and the total trajectory in the sample of an alpha particle having energy E=E_(max).

For a better explanation, according to FIG. II.2, there are two cases: E_(ind) less than the intermediate energy denoted (E_(m)) and E_(ind) greater than the intermediate energy E_(m).

1—case where E_(ind) is less than E_(m):

FIG. II-3 in the drawings section shows the variation in the sample space, from (θ_(i)) and (θ_(c)) as a function of the residual energy and trajectory (X), of an alpha particle of energy E_(i) and E_(ind)∈[E_(min), E_(m)] emitted of a radioactive core located at the height h from the surface of the detector. It is observed that:

-   -   when the residual energy decreases, (θ_(c)) decreases and         (θ_(i)) increases.     -   for (θ_(i))=0 (normal incidence for each position h), the         corresponding (θ_(c)) has the larger value for each position h,         and when (θ_(i)) increases, (θ_(c)) decreases. Thus, (θ_(ic))         for each position h is none other than the intersection between         the angle of incidence of the particle to be detected (θ_(i))         and the resolution angle of the detector (θ_(c)), because from         this value, (θ_(i)) is greater than (θ_(c)).

2—case where E_(ind) is greater than E_(m);

At first, in the interval of residual energy [E_(m), E_(ind)], it is observed that when the residual energy decreases, (θ_(i)) increases and (θ_(c)) increases.

Secondly, in the interval of residual energy [θ, E_(m)], it is observed that when the residual energy decreases, (θ_(i)) increases and (θ_(c)) decreases.

For (θ_(i))=0 (normal incidence for each position h), the corresponding angle (θc) is also greater than zero for each position h and when (θ_(i)) increases, (θ_(c)) increases and then decreases, thus when the function appears as (θ_(c))=F(E_(res)) {strong phase of decrease of (θ_(c)) for E_(res)∈[E_(m), E_(i)]}, (θ_(i)) will have one and only one intersection with (θ_(c)) for each position h, i.e. in the first interval or in the second interval of energy, from this value of intersection, (θ_(i)) is greater than (θ_(c)).

Thus, in a general manner, when the function appears as (θ_(c))=F(E_(res)), for each position h where E_(min)≦E_(ind)≦E_(max), (θ_(ic)) is none other than the intersection between the angle of incidence of the particle to be detected (θ_(i)) and the resolution angle of the detector (θ_(c)).

A-3—Case of Nuclei Emitters of Alpha Particles Having Energy E_(i) and for which E_(ind)>E_(max),

These are nuclei located in the volume for contribution to the tracks defined by h∈[0, R_(i)−R_(max)], for the case without an absorbent medium, where R_(max) and the total trajectory in the sample of an alpha particle having energy E=E_(max).

In the same manner as above, when the function has the appearance (θ_(c))=F(E_(res)), there will be either zero intersections, or 2 intersections between (θ_(i)) and (θ_(c)) for each (h) which will be denoted (θ_(ic1)) and (θ_(ic2)) which correspond respectively to the maximum energy and to the one that is the minimum detectable at oblique incidence for each (h).

Thus, for E_(ind)>E_(max), only the particles emitted at angles between (θ_(ic1)) and (θ_(ic2)) will be detected if they reach the detector.

B) Evaluation of the Function (θ_(ic)).

In order to determine (θ_(ic)), (θ_(ic1)), and (θ_(ic2)) as a function of h, for the case without an absorbent medium, proceed as follows: (θ_(ic))=f ₁(E _(res))=f ₂(x), E _(res) ∈E[0, E _(i) ], x∈E[0,R _(i)]

-   -   R_(i) is the total trajectory of the particle of energy E_(i) in         the medium analyzed.     -   x: is the trajectory of the particle of energy E_(i) in the         medium analyzed, this particle reaches the detector with an         energy E_(res), x will be determined from the relation between         energy and particle in the medium analyzed [2].         (θ_(i))=arccos(h/x)=g ₂(h,x)=g ₁(h,E _(res)), x∈[h,R _(i)]     -   h is the position of the radioactive nucleus emitter of the         particle of energy E_(i) relative to the surface of the         detector.     -   x: is the trajectory of the particle of energy E_(i) in the         medium analyzed and varies from h to R_(i) for each position h.

In fact, (θ_(ic)) will be determined by solving the equation: g₂(h,x)=f₂(x) for the heights h_(i) where it is confirmed that E_(min)≦E_(ind)≦E_(max), (θ_(ic1)) and (θ_(ic2)) will be determined by solving the equation: g₂(h,x)=f₂(x) for the heights h_(i) where it is confirmed that E_(ind)>E_(max). Thus, a determination is done of (θ_(ic))=f(h) for the values of h for which E_(min)≦E_(ind)≦E_(max), (θ_(ic1))=f_(a)(h) and (θ_(ic2))=f_(n)(h) for the values of h for which E_(ind)>E_(max).

It is possible to convert from a function f₁(E_(res)) to its equivalent f₂(x) for each type of alpha particle having energy (E_(i)) and total trajectory in the sample analyzed (R_(i)), emitted from a nucleus located in a position (h) within the sample to be analyzed, either by the slope of the function of Linear Energy Transfer (L.E.T.=dE/dx) of the alpha particles in a medium, or by the Energy-Trajectory relation [E=f(R)] of the alpha particles in a medium:

-   -   For the case without an absorbing medium:     -   Using the function of Linear Energy Transfer (dE/dx) in the         medium to be analyzed:         $E_{res} = {{\int_{0}^{R_{i} - x}\quad{\frac{\mathbb{d}E}{\mathbb{d}x} \cdot {\mathbb{d}x}}} = {E_{i} - {\int_{R_{i} - x}^{R_{i}}{\frac{\mathbb{d}E}{\mathbb{d}x} \cdot \quad{\mathbb{d}x}}}}}$         thus an E_(res) corresponds for each x and, as a result, the         conversion of a function f₁(E_(res)) to its equivalent f₂(x) and         vice-versa.     -   Using the Energy-Trajectory relation, [E=f(R)] in the medium to         be analyzed:

The function E=f(R) is used and the change in the variable (R) is made by (R_(i)−x), this function becomes: E_(res)=31 _(i)−x), thus for each x, the corresponding E_(res) will be determined and as a result, the conversion from one function f₁(E_(res)) to its equivalent f₂(x) and vice-versa.

-   -   For the case with an absorbing medium:     -   Using the Energy-Trajectory relation [E=f(R)]:

According to FIG. II-4 shown in the drawings section, where [(M1): medium to be analyzed. (M2): absorbent medium. (h_(a)): thickness of the absorbent medium. x₁ and x₂ are, respectively, the trajectory of the alpha particle in the analyzed medium and the one in the absorbent medium.]:

Given: ${{\cos\left( \theta_{i} \right)} = {\frac{h - h_{a}}{x_{1}} = {\frac{h_{a}}{x_{2}} = \frac{h}{x}}}},{x = {x_{1} + x_{2}}}$ i.e.:

-   -   E=f(R), the Energy-Trajectory function of alpha particles in the         medium to be analyzed.     -   R=g(E), the Energy-Trajectory function of alpha particles in the         absorbent medium.     -   E_(res1) and E_(res) are, respectively, the residual energy of         the incident alpha particle at the input of the absorbant medium         and that at the input of the SSNTD.

Given E_(res1)=f(R_(i)−x_(i)), thus R_(lab)=g(E_(res1))=g∘f(R_(i)−x₁). Where R_(lab) is the total trajectory of the alpha particle having energy E_(res1) in the absorbent medium. As a result, [R_(lab)−x₂]=[g∘f(R_(i)−x₁)−x₂]=g(E_(res)), or ${x_{2} = {{x_{1} \cdot \frac{h_{a}}{h - h_{a}}} = {x \cdot \frac{h_{a}}{h}}}},$ In fact, ${{g \cdot {f\left( {R_{i} - {x \cdot \frac{h_{a}}{h - h_{a}}}} \right)}} - {x \cdot \frac{h_{a}}{h}}} = {g\left( E_{res} \right)}$

As a result, for each (h), each x corresponds to a residual energy E_(res) at the input of the SSNTD. Thus, for each position (h_(i)), the function g₂(h_(i),x) will be transformed into its equivalent function g₁(h_(i),E_(res)) and vice-versa. As a result, the resolution for each (h_(i)) of the equation g₁(h_(i),E_(res))=f₁(E_(res)) will make it possible to determine θ_(ic)(h_(i)) for each (h_(i))∈[h_(o),R_(i)+h_(o)−f⁻¹∘g⁻¹(h_(o))].

-   -   Using the function for Linear Energy Transfer (dE/dx);

The calculation methodology is the same as the one using the Energy-Trajectory relation, only instead of the functions Energy-Trajectory and Trajectory-Energy in the analyzed medium and the absorbent medium, they are replaced with the concept of the integral of the function for Linear Energy Transfer as described in the case without absorbent medium.

II-1-2—Determination of the Probability P_(i)(r, h).

II-1-2-a—Case of Nucleus Emitting Alpha Particles of Energy E_(i) and for which E_(min)≦E_(ind)≦E_(max).

According to FIG. (II-5) shown in the drawings section, the proportion of particles emitted from a nucleus located in a position of coordinates (r, h) and at angles of incidence less than the critical angle of incidence (θ_(ic)) which reach the detector is the ratio of the surface of intersection {between the surface described by the angle (θ_(ic)) on the plane of the detector denoted s(θ_(ic)) and that of the detector denoted S(R)} and the surface described by the angle (θ_(ic)) on the plane of the detector denoted s(θ_(ic)). This proportion is described by the following equation: $\left\{ \frac{{S\left( \theta_{ic} \right)}\bigcap{S(R)}}{S\left( \theta_{ic} \right)} \right\}$

Thus, the probability that a particle of energy E_(i) emitted by the nucleus having index i contained in the volume dv and confirmed by E_(min)≦E_(ind)≦E_(max) can be detected by the solid state nuclear track detector is: ${{Pi}\left( {r,h} \right)} = {\left\{ \frac{1 - {\cos\quad\theta_{ic}}}{2} \right\}\left\{ \frac{{S\left( \theta_{ic} \right)}\bigcap{S(R)}}{S\left( \theta_{ic} \right)} \right\}}$

-   -   Evaluation of         $\left\{ \frac{{S\left( \theta_{ic} \right)}\bigcap{S(R)}}{S\left( \theta_{ic} \right)} \right\},$         with S(θ_(ic))=πR² _(c), where R_(c)=h·tg(θ_(ic)). (See FIG.         II-5)         a) for (θ_(ic))≧(θ₂), S(θ_(ic))∩S(R)=πR², where R: radius of         SSNTD         b) for (θ_(ic))≦(θ₂): 2 cases are present, (θ_(ic))≧(θ₁) and         (θ_(ic))≦(θ₁):     -   if (θ_(ic))≧(θ₁)         ${{S\left( \theta_{ic} \right)}\bigcap{S(R)}} = {2\left\lbrack {{\int_{r - {Rc}}^{I}{\sqrt{R_{c}^{2} - \left( {x - r} \right)^{2}}.\quad{\mathbb{d}x}}} + {\int_{I}^{R}{\sqrt{R^{2} - x^{2}}.\quad{\mathbb{d}x}}}} \right\rbrack}$         where $I = \frac{R^{2} + r^{2} - R_{c}^{2}}{2r}$     -   if (θ_(ic))≦(θ₁): S(θ_(ic))∩S(R)=S(θ_(ic))         II-1-2-b—Case of Nuclei Emitting Alpha Particles of Energy E_(i)         and for which E_(ind)>E_(max)

According to FIG. (II-5) shown in the drawings section, the proportion of particles emitted from a nucleus located in a position of coordinates (r, h) and at angles of incidence between (θ_(ic1)) and (θ_(ic2)) which reach the detector is the ratio of the surface of intersection {between the surface of the ring described by the angles (θ_(ic1)) and (θ_(ic2)) on the plane of the detector denoted s(θ_(ic1), θ_(ic2)) and that of the detector denoted S(R)} and the surface of the ring described by the angles (θ_(ic1)) and (θ_(ic2)) on the plane of the detector denoted s(θ_(ic1), θ_(ic2)). This proportion is described by the following equation: $\left\{ \frac{{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)}\bigcap{S(R)}}{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)} \right\}$

Thus, the probability that a particle of energy E_(i) emitted by the nucleus having index i contained in the volume dv and confirmed by E_(ind)>E_(max) can be detected by the solid state nuclear track detector is: ${{Pi}\left( {r,h} \right)} = {\left\{ \frac{{\cos\quad\theta_{{ic}\quad 1}} - {\cos\quad\theta_{{ic}\quad 2}}}{2} \right\}\left\{ \frac{{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)}\bigcap{S(R)}}{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)} \right\}}$ $◆\quad{Evaluation}\quad{of}\quad\left\{ \frac{{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)}\bigcap{S(R)}}{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)} \right\}$ with: −S(θ_(ic1),θ_(ic2))=π(R_(c2) ²−R_(c2) ²−R_(c1) ²) where R_(c1)=h·tg(θ_(ic1)), R_(c2)=h·tg(θ_(ic2))

-   -   −S(θ_(ic1),θ_(ic2))∩S(R)=S(θ_(ic2))∩S(R)−S(θ_(ic1))∩S(R)         a) for (θ_(ic1))≧(θ₂),S(θ_(ic1),θ_(ic2))∩S(R)=πR²−πR²=0     -   where R: radius of the SSNTD         b) for (θ₁)≦(θ_(ic1))≦(θ₂): there are 2 cases, (θ_(ic2))≧(θ₂)         and (θ_(ic2))≦(θ₂)     -   if (θ_(ic2))≧(θ₂):         ${{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)}\bigcap{S(R)}} = {{\pi\quad R^{2}} - {2\left\lbrack {{\int_{r - R_{c\quad 1}}^{I}{\sqrt{R_{c\quad 1}^{2} - \left( {x - r} \right)^{2}} \cdot \quad{\mathbb{d}x}}} + {\int_{I_{1}}^{R}{\sqrt{R^{2} - x^{2}} \cdot \quad{\mathbb{d}x}}}} \right\rbrack}}$         ${{where}\quad I_{1}} = \frac{R^{2} + r^{2} - R_{c\quad 1}^{2}}{2r}$         ◆if  (θ_(ic  2)) ≤ (θ₂):S(θ_(ic  1), θ_(ic  2))⋂S(R) = S(θ_(ic  2))⋂S(R) − S(θ_(ic  1))⋂S(R)         where:         ${{S\left( \theta_{{ic}\quad 2} \right)}\bigcap{S(R)}} = {2\left\lbrack {{\int_{r - R_{c\quad 2}}^{I_{2}}{\sqrt{R_{c\quad 2}^{2} - \left( {x - r} \right)^{2}} \cdot \quad{\mathbb{d}x}}} + {\int_{I_{2}}^{R}{\sqrt{R^{2} - x^{2}} \cdot \quad{\mathbb{d}x}}}} \right\rbrack}$         I₂ = R² + r² − R_(c  2)²  [sic] and         ${{S\left( \theta_{{ic}\quad 1} \right)}\bigcap{S(R)}} = {2\left\lbrack {{\int_{r - R_{c\quad 1}}^{I_{1}}{\sqrt{R_{c\quad 1}^{2} - \left( {x - r} \right)^{2}} \cdot \quad{\mathbb{d}x}}} + {\int_{I_{1}}^{R}{\sqrt{R^{2} - x^{2}} \cdot \quad{\mathbb{d}x}}}} \right\rbrack}$         $I_{1} = \frac{R^{2} + r^{2} - R_{c\quad 1}^{2}}{2r}$         c) for (θ_(ic1))≦(θ_(i)): there are 2 cases, (θ_(ic2))≧(θ₁) and         (θ_(ic2))≦(θ₁):         1—for (θ_(ic2))≧(θ₁): there are 2 cases, (θ_(ic2))≧(θ₂) and         (θ_(ic2))≦(θ₂):     -   if (θ_(ic2))≧(θ₂): S(θ_(ic1), θ_(ic2))∩S(R)=πR²−πR_(c1) ²         ${{{◆if}\quad\left( \theta_{{ic}\quad 2} \right)} \leq {{\left( \theta_{2} \right)\text{:}{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)}}\bigcap{S(R)}}} = {{2\left\lbrack {{\int_{r - R_{c\quad 2}}^{I_{2}}{\sqrt{R_{c\quad 2}^{2} - \left( {x - r} \right)^{2}} \cdot \quad{\mathbb{d}x}}} + {\int_{I_{2}}^{R}{\sqrt{R^{2} - x^{2}} \cdot \quad{\mathbb{d}x}}}} \right\rbrack} - {\pi R}_{c\quad 1}^{2}}$         2-  for  (θ_(ic  2)) ≤ (θ₁):S(θ_(ic  1), θ_(ic  2))⋂S(R) = S(θ_(ic  1), θ_(ic  2)) = π  R_(c  2)² − π  R_(c  1)²           II-2—Generalization for the Total Contribution Volume for         Resolution of Tracks:

According to the above, the total number of tracks observed and which are due to the radioactive nuclei of index i contained in the total volume: n _(it)=2π·A _(ci) ^(te) ·∫∫P _(i)(r,h)·rdrdh

For the case E_(min)≦E_(ind)≦E_(max), a comparison of (θ_(ic)) with (θ₁) and (θ₂) allows to determine respectively r₁(h) and r₂(h) where these latter are characteristics of respectively the radii defining the interval of the volume where (θ_(ic))≦(θ₁) from the one where (θ_(ic))≧(θ₁) and the ones defining the interval of the volume where (θ_(ic))≦(θ₂) from the one where (θ_(ic))≧(θ₂).

For the case E_(ind)>E_(max), a comparison of (θ_(ic1)) with (θ₁) and (θ₂) allows to determine respectively r_(1,1)(h) and r_(2,1)(h) where these latter are characteristics of respectively the radii defining the interval of the volume where (θ_(ic))≦(θ₁) [illegible] (θ_(ic))≧(θ₁) and the ones defining the interval of the volume where (θ_(ic1))≦(θ₂) [illegible] (θ_(ic1))≧(θ₂). And a comparison of (θ_(ic2)) with (θ₁) and (θ₂) allows us to determine respectively r_(1,2)(h) and r_(2,2)(h) where these latter are characteristics of respectively the radii defining the interval of the volume where (θ_(ic2))≦(θ₁) from the one where (θ_(ic2))≧(θ₁) and the ones defining the interval of the volume where (θ_(ic2))≦(θ₂) from the one where (θ_(ic2))≧(θ₂).

The next part of the procedure is an integration by parts of the interval over r and h.

To determine r₁(h), r₂(h), r_(1,1)(h), r_(2,1)(h), r_(1,2)(h), and r_(2,2)(h), it is simply a matter of solving the following equations: It has: −(θ_(ic))=f(h); (θ_(ic1))=f_(a)(h); (θ_(ic2))=f_(b)(h) ${{- \left( \theta_{1} \right)} = {{arc}\quad{\cos\quad\left\lbrack \frac{h}{\sqrt{h^{2} + \left( {R - r} \right)^{2}}} \right\rbrack}}};{\left( \theta_{2} \right) = {{arc}\quad{\cos\left\lbrack \frac{h}{\sqrt{h^{2} + \left( {R + r} \right)^{2}}} \right\rbrack}}}$ thus:

-   -   *(θ_(ic))=(θ₁)         r₁(h); *(θ_(ic))=(θ₂)         r₂(h)     -   *(θ_(ic1))=(θ₁)         r_(1,1)(h); *(θ_(ic1))=(θ₂)         r_(2,1)(h)     -   *(θ_(ic2))=(θ₁)         r_(1,2)(h); *(θ_(ic2))=(θ₂)         r_(2,2)(h)

Actually, the total volume of the sample contributing to the tracks (V_(t)) will be the total of a series of the volume V_(f) (V_(t)=ΣV_(j)), where each V_(j) itself corresponds to a probability P_(i,j)(r,h).

In conclusion and in a general manner, the total number of tracks observed and which are due to the radioactive nuclei of index i contained in the total volume: $n_{it} = {A_{ci}^{te} \cdot {\sum\limits_{j}{\int{\int{\int{{P_{ij}\left( {r,h} \right)} \cdot {\mathbb{d}v_{j}}}}}}}}$ where: Pij(r,h) and vj will be defined as a function of r₁(h), r₂(h), r_(1,1)(h), r_(2,1)(h), r_(1,2)(h), and r_(2,2)(h) with: ${{{Pi}\left( {r,h} \right)} = {\left( \frac{1 - {\cos\quad\theta_{ic}}}{2} \right)\left\{ \frac{{S\left( \theta_{ic} \right)}\bigcap{S(R)}}{S\left( \theta_{ic} \right)} \right\}}};{{{for}\quad E\quad\min} \leq {E\quad{ind}} \leq {E\quad\max}}$ ${{{Pi}\left( {r,h} \right)} = {\left( \frac{{\cos\quad\theta_{{ic}\quad 1}} - {\cos\quad\theta_{{ic}\quad 2}}}{2} \right)\left\{ \frac{{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)}\bigcap{S(R)}}{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)} \right\}}};{{{for}\quad E\quad{ind}} > {E\quad\max}}$ where: $\left\{ \frac{{S\left( \theta_{ic} \right)}\bigcap{S(R)}}{S\left( \theta_{ic} \right)} \right\}\quad{and}\quad\left\{ \frac{{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)}\bigcap{S(R)}}{S\left( {\theta_{{ic}\quad 1},\theta_{{ic}\quad 2}} \right)} \right\}$ are well defined for c [sic] III—Application Example:

Taking the simplest case, involving the media where the total trajectory of the particle α_(i) of energy E_(i) is much less than the radius of the detector (R_(i)<<R) {case most reported due to the fact that this condition is verifiable for alpha particles in almost all liquid and solid samples} and a detector having a response function for conditions of developing and reading the well-defined tracks such as E_(min)=0 and E_(max)>E_(j).

Following the above gives: ${\cos\quad\left( \theta_{2} \right)} = \frac{h}{\sqrt{h^{2} + \left( {R + r} \right)^{2}}}$ Thus, (θ_(2min)) for (Ri<<R) approaches $\frac{\pi}{2},$ and as a result, (θ_(ic)) is always less than (θ₂)

Actually, a comparison of (θ_(ic)) with (θ₁) allows us to determine r₁(h) and, as a result: $n_{it} = {2\quad{\pi \cdot A_{ci}^{te} \cdot \left\{ {{\int_{0}^{R_{1^{\prime}i}}{\int_{0}^{(h)}{\frac{1 - {\cos\quad\theta_{ic}}}{2}r\quad{\mathbb{d}r}\quad{\mathbb{d}h}}}} + {\int_{0}^{R_{1}}{\int_{r_{1}{(h)}}^{R}{\frac{1 - {\cos\quad\theta_{ic}}}{2}\left( \frac{{S\left( \theta_{ic} \right)}\bigcap{S(R)}}{S\left( \theta_{ic} \right)} \right)r\quad{\mathbb{d}r}\quad{\mathbb{d}h}}}}} \right\}}}$ where: _S(θ_(ic)) = π  R_(c)², R_(c) = h ⋅ tg(θ_(ic)) ${{{\_ S}\left( \theta_{ic} \right)}\bigcap{S(R)}} = {2\left\lbrack {{\int_{r - R_{c\quad}}^{I}{\sqrt{R_{c\quad}^{2} - \left( {x - r} \right)^{2}} \cdot \quad{\mathbb{d}x}}} + {\int_{I}^{R}{\sqrt{R^{2} - x^{2}} \cdot \quad{\mathbb{d}x}}}} \right\rbrack}$ $I = \frac{R^{2} + r^{2} - R_{c}^{2}}{2r}$ IV—Solution to the Problem Involving Media Having Several Types of Nuclei Emitting Particles with Different Energies and Activities.

A technique for extracting energy is used whereby an absorbent medium must be installed in contact with the solid state nuclear track detector in a manner such that only those particles having the largest energy can reach the detector and the other energies follow in sequence.

And if it appears that for a given medium, the ability to separate the energies by the SSNTD using an absorbent medium is not feasible, then an absorbent medium of thickness (h_(a)) is installed in contact with the detector with different thicknesses (h_(a)) of the absorbent medium, in a manner such that the number of cases of thicknesses (h_(a)) of the absorbent medium coincides with the total number of different alpha emitters of the medium to be analyzed. As a result, one has the number of equations linking the total number of tracks observed for each (h_(a)) to the activities of the different alpha emitters, [illegible] of the alpha activities of the different emitters (the unknowns). Thus, [illegible] alpha activities of the different emitters.

For the determination of (θ_(ic)(h)), (θ_(ic1)(h)), (θ_(ic2)(h)) of each type of nucleus in the presence of an absorbent medium, the method is to use two equations (θ_(c)=f(E_(res))) and (θ_(i)=f(E_(res))) (see part B of paragraph II-1-1) and the calculation methodology is the same while taking into account the presence of the absorbing medium.

V—Conclusion:

This direct technique allowing the dosimetry of the radioactive nucleii in different types of samples has been developed at home and in an individual manner. It is based on the use of a novel method adapted to the experimental conditions using the solid state nuclear track detectors.

It consists in calculating the probability P_(i)(r,h) that an alpha particle emitted by a radioactive nucleus contained in the volume dv of a sample can be detected by a solid state nuclear track detector, where (dv) is the volume element of a cylinder having a radius between r and r+dr and a height between h and h+dh inside the sample to be analyzed.

This method is easy to implement even though its formalism appears quite complex. It is precise, absolute, and not costly.

Solely for reasons of simplification of the calculations the radius is set of the detector equal to the radius of the irradiation cell, the methodology of calculation is the same for a radius of the detector less than the radius of the irradiation cell with small modifications at the step of evaluation of the surfaces of intersection {S(θ_(ic))∩S(R)} and {S(θ_(ic1),θ_(ic2))∩S(R)}, this time, the calculation must be made of the surfaces of intersection for r (characterizing the position of the emitting nucleus) less than the radius of the detector and for r greater than the radius of the detector.

In a general manner, regardless of the appearance of the experimental function of the intrinsic critical angle of a detector θ_(c)=F(E_(res)) and for the position (h_(i)), more than 2 intersections with θ_(i)(h_(i)) are present, the efficiency of the detection for these positions (h_(i)) will be a series of combinations between the two general formulas for efficiencies of detection described in this paper, the one corresponding to a single intersection and the other corresponding to two intersections, where the nature of this combination depends on the number of intersections.

REFERENCES

-   [1] G. Somogyi, Nucl. Instr. Meth., 173, 21 (1980). -   [2] J. F. Ziegler, Programme TRIM. 

1. A method for dosimetry of radioactive nuclei in different types of samples, comprising the steps of: using solid state nuclear track detectors wherein implementation allows determination of a function defining actual sensitivity of detection of a solid state nuclear track detector for each type of nucleus emitting alpha particles having energy E_(i) located in each position within a sample to be analyzed, being based on a actual function of critical angles of incidence (θ_(ic)), (θ_(ic1)) and (θ_(ic2)).
 2. The method for dosimetry of claim 1, further comprising: identifying (θ_(ic)), (θ_(ic1)) and (θ_(ic2)) by demonstration of relation between the angle of incidence of the particle to be detected (θ_(i)) and the resolution angle of the detector (θ_(c)).
 3. The method of dosimetry of claim 2, further comprising: evaluating (θ_(ic)), (θ_(ic1)) and (θ_(ic2)) wherein (θ_(ic))=f(h) for the h such that E_(min)≦E_(ind)≦E_(max), (θ_(ic1))=f_(a)(h) and (θ_(ic2))=f_(b)(h) for the h such that E_(ind)>E_(max).
 4. The method for dosimetry of claim 3, further comprising: identifying and evaluating surfaces of intersections {S(θ_(ic))∩S(R)} and {S(θ_(ic1), θ_(ic2))∩S(R)} representing respectively the surface of intersection between the surface described by the angle (θ_(ic)) on the plane of the detector and the surface of the detector of radius (R), and the surface of intersection between the surface described by the ring defined by (θ_(ic1)) and (θ_(ic2)) on the plane of the detector and the surface of the detector of radius (R), wherein outcome of the determination of the function representing the efficiency of the detection P_(i)(r,h).
 5. The method for dosimetry of claim 4, further comprising: generalizing to total volume by the separation of the total volume of contribution to the tracks (V_(t)) into a series of volume V_(j)(V_(i)=ΣV_(j)), where each V_(j) itself corresponds to a function representing the efficiency of detection P_(i,j)(r,h). 